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ABSTRACT 



In this investigation, use was made of an electronic analog 
for simulating the operation of a free-piston gas generator* The 
basic concept of the analog was developed in a previous investiga- 
tion and involved the creating of the various free-piston forces as 
functions of displacement and integrating these combined forces in 
accordance with Newton’s second law of motion to obtain a solution 
representing piston motion - all by electronic analog means. 

Modifications were made to the basic analog circuit of the pre- 
vious investigation to overcome an unstable operating condition. 

The modified circuit was then employed for determining performance 
characteristics of an actual free-piston gas generator under varied 
load conditions. 

The results of these studies showed that the operation of a 
free-piston gas generator could be simulated to a reasonable degree 
of accuracy and stability by use of an electronic analog and that 
this method might be a valuable means for design and performance 
predictions of free-piston engine systems. 
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NOMENCLATURE 



C 

C 

f r 

F fr 

h 

I DP 
k 

LHV 

m d 

m 



rn f 



M 



- Potentiometer constant, dimensionless 

2 



English Letter Symbols 
a 
A 

b 

BSFC 



- Cross-sectional area, ft 

- Bounce cylinder clearance with engine pistons in contact, ft 

- Brake specific fuel consumption, Ibm fuel/shphr 



n’ 



- Compressor cylinder clearance with engine pistons in 
contact, ft 

- Coulomb friction force, lbf 

- Electrical capacitance, fd 

- Frequency, cycl^s/min 

- Frictional force, lbf 

- Enthalpy, BTU/lbm 

- Inner dead point piston position, ft 

- Ratio of specific heats, dimensionless 

- Lower heating value of fuel, BTU/lbm 

- Mass of air delivered by the compressor, Ibm 

- Mass of engine intake air, Ibm 

- Mass of fuel, Ibm 

- Mass of piston, Ibm 

- Polytropic exponent for bounce cylinder compression and 
expansion processes, dimensionless 

- Polytropic exponent for compression cylinder compression 
and expansion processes, dimensionless 

- Polytropic exponent for engine cylinder compression 
process, dimensionless 

- Polytropic exponent for engine cylinder expansion 
process, dimensionless 
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English Letter Symbols (continued) 



ODP 



P 

o 




V 

P 

Wk 



x 



X 

p 



- Outer dead point piston position, ft 

2 

- Ambient air pressure, lbf/ft 

2 

- Compressor intake pressure, lbf/ft 

2 

- Compressor discharge pressure, lbf/ft 

2 

- Scavenge air receiver pressure, lbf/ft 

2 

- Engine cylinder pressure during scavenging, lbf/ft 

2 

- Gas delivery pressure, lbf/ft 

2 

- Maximum or initial bounce cylinder pressure, lbf/ft 

- Energy of fuel^BTU 

- Gas constant for air, ft-lbf/lbm-°R 

- Electrical resistance, ohms 

- Stroke, ft 

- Shaft horsepower, hp 

- Time, sec 

- Ambient air temperature, °R 

- Compressor discharge temperature, °R 

- Engine exhaust gas temperature, °R 

- Engine intake air temperature, °R 

- Gas delivery temperature, °R 

3 

- Engine cylinder volume at point of port closure, ft 

- Work, ft- Ibf 

* Engine piston displacement from midpoint or point of 
contact of pistons, ft 

- Engine piston displacement at port closure, ft 
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English Letter Symbols (continued) 

>tj, y s, z - Piston velocity, ft/sec 

2 

x - Piston acceleration, ft/sec 

y - Compressor piston displacement from compressor cylinder 

head, ft 

- Length of air intake portion of compressor stroke, ft 
z - Bounce piston displacement from bounce cylinder he*d, ft 

Subscripts : 

b - Refers to bounce cylinder 

c - Refers to compressor cylinder 

e - Refers to engine cylinder 

f - Refers to operational amplifier feedback 

s - Refers to isentropic process 

Greek Letter Symbols: 

i - Magnitude scaling factor, units of i/volt 

t - Time scaling factor, computer time/real time 

A * Difference, dimensionless 

£ - Relay bias voltage, volts 

Vy - Efficiency, percent 

Miscellaneous : 

i - Scaled computer voltage quantity - related to real 

physical quantity, i, by scale factor, o' - thus 

i = ©<, I 
1 
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1. Introduction 



An electronic analog of the operation of a free-piston gas gener- 
ator or gasifier based on the "spring-mass" nature of the free-piston 
engine system was undertaken by LT Arthur E. PLOW, USN. OJ m that 
investigation, the first law of thermodynamics was solved on an analog 
computer three concurrent times representing electronic analogs of an 
internal combustion engine, a reciprocating compressor, and a bounce 
cylinder "gas spring". Friction forces were also introduced. The re- 
sults of these four analogs were combined and integrated in accordance 

with Newton* s second law of motion to yield a solution representing 

•% 

piston motion. In addition a comparison of analog computer solution 

results and actual performance data at approximate rated output condi- 

2 

tions was also made for the SIGMA organization of France model GS-34 
free-piston gasifier. 

The results of LT PLOW®s investigation showed that the various 
nonlinear "gas-spring" forces of the gasifier (engine, compressor and 
bounce cylinder forces) as well as a friction force could be simulated 
by electronic analog methods. The operation of the bounce cylinder 
analog, however, was non-stable or drifting and tended to disrupt the 
problem. The innaccuracy of this analog resulted from errors in a 
two-quadrant electronic multiplication performed in this circuit, 
these errors being due to the inherent inability of an electronic mult 
plier to give uniform results in multiplications involving more than 
one quadrant (i.e. both inputs having the same sign). 

^Numerals in brackets refer to bibliography, 

2 « 0 e t B » 

Societe Xndustrielle Generate de Mecanique Appliquee (SIGM/)» 
Vtnissieux s France. 
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Relatively few techniques for predicting the design and performance 
of free-piston engine systems may be found in the available literature. 
These, for the most part, involve rather complex numerical processes 
as illustrated by London’s [ 2 ~] thermodynamic-dynamic analysis of a 
free-piston gas generator system in which a calculated "basic design’* 
may be extrapolated by certain affinity relations to different plant 
capacities and dimensions. 

In this investigation the electronic analog method as developed by 
PLOW {Jl] but modified for stability of operation was employed for deter- 
mining performance characteristics of a free piston gas generator under 

* 

varied load conditions. To this end full and partial load studies were 
made of the SIGMA GS-34 gasifier by introducing appropriate fixed geo- 
metric and variable operating parameters to the analog problem. This 
procedure could be adapted as well to performance predictions of modi- 
fied gas generator designs and thus would be of considerable benefit in 
the development of different sized machines. 

In the available literature there is little information on predic- 
tion of part load behavior of a free-piston gas generator. The analog 
method might well be a useful means in this type of analysis. 
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2, Characteristics of the Free-Piston Engine System* 



The free-piston engine system consists of an opposed piston 9 super- 
charged two- stroke Diesel engine driving a reciprocating air-compressor 
and may be combined with a turbine which utilizes the expansion of the 
exhaust gases down to an atmospheric condition. The engine reciprocat- 
ing work is directly used to provide the compressor work requirement. 

The free piston system is characterized by (a) constructional simplicity 
due to the absence of cranks and bearings; (b) almost complete freedom 
from vibration provided by th crankless opposed piston design; and (c) 
the absence of piston-cylinder side thrust and therefore reduction of 
cylinder wear as introduced by crank-connecting rod systems. 

As shown in Fig. 1, there are two basic applications of the free- 
piston system. Fig. 1 (a) represents an internal-combustion engine air 
compressor combination (free-piston air compressor) in which the useful 
output is compressed air foi pneumatic purposes. Fig. 1 (b) shows an 
air compressor internal combustion-engine combination (free-piston gas 
generator or gasifier) for the production of hot gases under pressure. 

The gasifier illustrated has an inboard compression, outboard bounce, 
and common central combustion chamber. Useful shaft work is derived 
from the hot pressurized gas by expansion to atmospheric pressure 
through a turbine. This investigation involved the latter configuration, 
i.e. s the free piston gas generator or gasifier. 
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(a). FREE -PISTON COMPRESSOR SYSTEM. 




TO ATMOSPHERE 

(b). FREE -PISTON GAS GENERATOR-TURBINE SYSTEM. 

Figure 1. Basic Applications of the Free-Piston System. 
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Referring to Fig, 1 (b), the gasifier consists of two opposed 
pistons having equal and symmetric strokes. The Diesel or power cylin- 
der is in the center and operates as a two-stroke Diesel engine super- 
charged to a pressure of several atmospheres (according to turbine lead 
condition). The two single-acting compressor cylinders are located on 
both ends of the central housing or scavenge air receiver. The cushion 
or bounce cylinders which store the energy for the return stroke are 
located at the outboard ends of the gasifier. Fresh air is taken in 
at atmospheric pressure through the air intake valves and is discharged 
into the scavenge air receiver surrounding the power cylinder. This com- 
pressed air is used for scavenging and charging of the power cylinder. 

The hot combustion gases, mixed with the excess scavenging air, ex- 
haust into the receiver and produce the useful power in the gas turbine. 
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3. Dynamics of Piston Motion. 



In the free-piston engine, the absence of a crank results in a 
"spring-mass" system which operates at a natural frequency depending 
on the mass of the piston assemblies and the nature of the nonlinear 
"gas springs" within the cylinders. London's Q2J spring-mass analog 
for the free-piston gas generator configuration is shown diagrammatically 
in Fig. 2 below: 




Inlet Port Exhaust Port 



Figure 2. Spring-Mass Analog for Free-Piston Gas-Generator System. 

Pressure or force versus volume or displacement diagrams for the 
power, compressor and bounce cylinders corresponding to the nonlinear 
spring forces of Fig. 2 are shown in Fig. 3. The pressure volume dia- 
gram of the power cylinder is representative of the Diesel cycle which 
was employed in this investigation. Also shown in this figure is a 
constant friction force that is independent of piston speed and gas 
pressures. Justification for this assumption may be found in Bobrowsky 

li2. 

During the outward stroke of the pistons the energy produced in 
the power or engine cylinder in addition to that of re-expansion of 
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Power Cylinder 
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Compressor Cylinder 
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(c) Friction Force 



Volume or Displacement 
(d) Bounce Cylinder 



Figure 3, Forces Acting on the Free Piston, 
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air in the clearance space of the compressor cylinders less friction 



losses is stored in the cushion or bounce cylinders. 



Wlc 



engine 
exp stroke 



= Wk 

bounce cyl 
compr stroke 



Wk + Wk a 0 

compressor friction 

exp stroke per stroke 



.3.1 



During the inward stroke the energy stored in the bounce cylinders is 
returned to the system again discounting friction losses, compressing 
and discharging the air in the compressor cylinders and compressing the 



air used for combustion in the power cylinder. 



Wk = Wk + Wk + Wk 

bounce cyl engine compressor friction 

exp stroke compr stroke compr stroke per stroke 

On a net work per cycle basis, 



3.2 



Wk 0 * Wk + Wk r . . 3.3 

engine compressor friction . 

cycle cycle per cycle 

For each half cycle the component work terms must balance since the 

kinetic energy in the piston mass is zero at the inner and outer dead 

points where piston velocity is zero. 

Newton's Second Law relates piston resultant force F and the 

R 

piston rate of change of momentum where the resultant force F may be 

R 

evaluated as a function of piston position from 

F_ (out stroke) = F + F - F, - F, 

R eng compr bounce fr 

3.4 

F_ (in stroke) = F, - F - F F . 

R bounce compr eng fr 

The vjork done by the resultant force acting on the piston is equal to 

the gain of kinetic energy of the piston. 
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4 0 General Electronic Analog of the Free-Piston Gas Generator* 

In the previous section it was shown that the free-piston engine 
may be considered as a “spring-mass' 1 system with a single degree of 
freedom wherein the nonlinear “gas spring" forces are of a thermodynamic 
nature* A constant static friction force opposing piston motion is 
also present* As previously developed in Plov; D3 , the differential 
equation of motion of the piston in a gas generator is a mathematical 
expression of Nex^ton's Second Law (Force = mass x acceleration) and 
is given by 

F e (x,x) +’ F c (x,x) - F b (x,x) - F fr (x) = Mx, 4.1 

where M is the piston mass 

(x,x) is the engine cylinder force on the piston 

(x,x) is the compressor cylinder force on the piston 

F, (x,x) is the bounce cylinder force on the piston 
b 

F- (x) is the friction force 

fr 

x is the piston displacement 

x is the piston velocity 

x is the piston acceleration 

It will be seen in later sections that the differential equations from 
which the engine, compressor and bounce cylinder forces or pressures 
are obtained are functions of piston displacement and velocity and that 
the friction force is dependent on piston velocity* Diagrams illustrat- 
ing the r gas spring" and friction forces acting on the piston are shown 
in Fig* 3* 

The electronic analog computing arrangement (from Plow QlJ) for the 



solution of equation 4,1 is shown in Fig, 4, The circuitry of the 
blocks representing the engine, compressor, bounce and friction ana- 
logs will be described in detail in following sections. Useful re- 
sults from the computing arrangement of Fig. 4 for determination of 
gas generator performance data will consist of piston displacement as 
a function of time and pressure-displacement diagrams for the engine, 
compressor, and bounce cylinders. Pressures and displacements will 
appear in terms of scaled voltages. 
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Figure 4. Basic Electronic Analog Computer Arrangement for 
Free-Piston Gas Generator. 
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5* System Inputs* Parameters* and Performance Relations* 

Inputs to the system will consist of: 
lo Fixed geometric parameters 0 

2* Initial piston displacement position (at outer dead point) 
3« Various cylinder pressure conditions* which are: 
a c Compressor intake and delivery pressures 0 
b 0 Engine pressure at start of compression 0 
Co Initial bounce cylinder pressure 0 
4 0 Fuel quantity 0 

So Thermodynamic* parameters such as polytropic exponents* 

Sol Geometric parameters* 

Variables and fixed quantities depending upon the geometry of the 
free-piston system are as follows: 

1* Engine piston area* A^ 0 

2 0 Compressor piston area* A^* 

3* Bounce piston area* A^« 

4* Engine piston displacement from midpoint or point of 
contact of pistons* x* 

So Compressor piston displacement from compressor cylinder 
head* y* 

6* Bounce piston displacement from bounce cylinder head* z* 

7o Engine piston displacement at port closure* 

8* Piston mass* M* 

The following relations also hold for the above parameters: 
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A b - 


A + A 9 

e c 


5.1.1 




y s 


X + C 8 


5.1,2 




Z as 


b - x 9 


5. 1.3 


and 


dx « 


dy s - dz 


5.1.4 



dt dt dt 

The quantities c and b in the above equations represent the clear- 
ances of the compressor and bounce cylinders respectively when the 
pistons are in contact at the midpoint or center of the machine* Fig* 

5 below shows these relations for one-half of the gas generator configura- 
tion* 



— 1 


1 C 


T 

J2 

1 

_L 


— x — • 


1 

J 1 












» I 

t 

Inlet Port 





x y z 



Figure 5* Illustration of Displacement and Clearance 
Parameters for Pistons in Contact at Mid- 
point* 



5*2 Pressures* 



In the free-piston gas generator ambient air of about one at- 
mosphere is compressed in the compressor to several atmospheres abso- 
lute pressure depending on turbine load condition and is delivered to 
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the scavenge air receiver during the inward stroke of the pistons. 

During the latter part of the power or outward stroke and commencement 
of inward stroke of the pistons, when both the exhaust ports and the in- 
take or scavenging ports are open, the scavenging air flows through the 
scavenge ports and through the full length of the power cylinder, forc- 
ing the exhaust gases through the exhaust ports and into the turbine in- 
let receiver (See Fig* 1), Some of the scavenging air follows the ex- 
haust gases into the receiver for good scavenging, but a portion of 
this supercharge air remains within the power cylinder and is compress- 
ed for the next combustion cycle during the inward stroke of the piston. 
Intake of fresh air into the compressor cylinders occurs during the out- 
ward piston stroke. 

The significant pressures in the free piston system are thus: 

1. P - ambient air pressure, 

2. - compressor intake pressure, 

3. P, - compressor discharge pressure, 

d 

4. P scavenge air receiver pressure, 

s 

5. P engine cylinder pressure during scavenging, 

P 

6. P - operating pressure at turbine inlet. 

Valve pressure drop allowances of five percent of upstream abso- 
lute pressure were arbitrarily chosen for each port passage of the gas 
generator. These include the compressor intake and discharge valves and 
the engine scavenge and exhaust ports. The ratio of downstream to up- 
stream pressure (absolute) for the compressor valves and engine ports is 
therefore equal to 0,95, 
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5.3 Temperatures 

As will be seen in succeeding subsections* the significant tempera- 
tures in the free-piston problem for the determination of the ir.asses of 
engine intake air, m^, and compressor air, and for the determina- 

tion of generator gas delivery temperature, T , are as follows: 

1. Ambient air temperature, T^o 

2. Engine intake air temperature, T^, assumed equal to com- 
pressor discharge temperature, T^. 

3. Engine exhaust gas temperature, T . 

5.4 Masses of Engine Intake Air and Compressor Discharge Air. 

Gas generator performance characteristics require the determination 
of engine intake and compressor discharge air masses. 

Employing the equation of state of a perfect gas, the mass of 
engine intake air m^, will be 



Pp Vp 



R T p 



5.4.1 



where is the engine cylinder volume at the point of port closure 

and R the gas constant for air (R = 53.3 ft-lbf/lbm - °R). 

Assuming engine intake air temperature, T , equal to compressor dis- 
charge temperature, T^, and considering compression and expansion pro- 
cesses in the various cylinders to be polytropic then 



rr\z = 




5.4.2 
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where is the polytropic compression exponent in the compressor 0 

Again employing the Equation of state of a perfect gas 9 the mass 

of air delivered by the compressor per cylinder per cycle 9 is given 

d 

by 

P. A Y. 

m , - 1 c 1 5.4.4 

d Tt ' 

o 

where is the area of the compressor cylinder and Y^ is the length 

of the air intake portion of the compressor stroke. 

The above expressions for engine intake and compressor discharge 
air masses involve only one -half of the free-piston gas generator and 
would have to be doubled to obtain total masses for both cylinders or 
sides of the system. 

5.5 Power Output of the System. 

In the actual free-piston gasifier-turbine combination hot combus- 
tion gases mixed with excess scavenging air exhaust from the gas genera- 
tor into the turbine receiver and then flow through the gas turbine which 
alone produces useful power. The thermodynamic state of the gases at the 
turbine inlet is dependent on the pressure and temperature of the exl aust 
gases. The product of the isentroptc available energy and the gas flow 
rate results in the rate at which work can be obtained from the exhaust 
gases. By multiplying this quantity by the turbine isentropic efficiency 
the effective power output at the turbine shaft finally can be determined. 

The exhaust gas temperature of a reciprocating internal combus- 
tion engine system is substantially less than the working substance tem- 
perature at the end of the power stroke. London considers the engine 
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exhaust ~ s g tempt Saturn eqrj dl r o the » f*« rxture of the 

r biow-dowr T period and has derive the to Howl.- i. .* engine e <: 

haust gas ter putafcure equation for a constant c "linden v. exhaust* 



e - jr [ L + ( k - 1) ^J ■ 5.5.1 



Tf = 



Subscripts 5 and 6 represent the state points at beginning and end of 
exhaust or "blow-down" period (See Fig 0 3), Temperature at state 5 is 
determined from the equation of state of a perfect gas* 

T* = g£ 

' s m e Q ; 5o5 ° 2 

where is mass of engine intake air as obtained from equation 

5o4„l and R is the gas constant for air. 

The generator gas delivery temperature* T „ is then equal to 
the mass average between engine exhaust gas temperature, and 

scavenge gas or engine intake air temperature* and is given by* 



T t = Te + ( 1 






5.5.3 



where m /m is the ratio of engine air eo compressed air delivered, 
e c 



m 



e 



and 



m, being obtained from equations 5 0 4 0 1 and 5o6 0 4 respectively* 



and assumed equal to compressor discharge temperature* is obtain- 

ed from 



T p - T d ~ To 



Pd 

Pi 



n c - 1 



n c 



5.5.4 



with quantities as defined in preceding subsections, 
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With the working substance considered a perfect gas, due to re- 
latively low delivery pressures of the free piston gasifier, enthalpy 
of exhaust gases may be considered a function of temperature alone 
and isentropic available energy determined,, The power out of the gas 
turbine is the product of the isentropic available energy, total exhaust 
gas flow rate considering the two sides of the gas-generator, and tur- 
bine isentropic efficiency,, Including appropriate conversion factors, 
the equation representing power output in horsepower is given by. 



( s hp) tur bine 



/a ^ S X 1 1 )L tLkJ 5.5.5 

33,000 ? 



whete^ h 



m , 



and 



is the isentropic available energy (Btu/lbm), 

is the mass (Ibm) of air delivered by the compressor per 

cylinder per cycle, 

is the piston frequency in cycles per minute, 
is the turbine isentropic efficiency. 
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6 0 Electronic Analog of the Two-Stroke Standard Diesel Cycle 0 

Fig* 3(a) shows a representative pressure versus? volume indicator 
diagram for the two-stroke air standard Diesel cycle assumed in this 
invest ig&t ion 0 With the absence of a crank mechanism in the free- 
piston gas generator* the stroke of the system and hence the engine com- 
pression ratio and compressor and bounce cylinder clearance volumes are 
all variable quantities dependent on the load conditions The only fix- 
ed geometrical point for the piston in the engine cylinder is that of 
port closureo The thermodynamic state of the engine intake air inass at 
this point can be determined for a given compressor discharge pressure 
which is also dependent on load conditions,, 

As previously considered, the working substances in all the cylinders 
of the free-piston system are assumed to be perfect gases and compression 
and expansion processes polytropic 0 Again referring to the indicator dia- 
gram of Figo 3(a), the two-stroke standard Diesel cycle can be seen to 
consist of the following processes* 

lo Polytropic compression* at exponent n of engine charge 
after port closure* 

2 0 Injection and burning of fuel to give addition of energy 
at constant pressure during first portion of power stroke* 

3o Polytropic expansion, at exponent n ^ 9 of engine charge and 
products of combustion* 

4* Engine "blow- down"' of exhaust gases at constant volume* 

5o Scavenging at constant pressure,. 

For a polytropic process involving a perfect gas 

PV° * constant, 6*1 
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where a is the polytropic exponent, by differentiation and con- 



sidering time as the independent variable, we obtain 

JP 



d P 

With tne volume given by 



^ d v 



V 



V — A e X ) 



and 



d V a d K , 

= Az d* = Ae “’ 



6.2 



6.3 



6.4 



6.5 



d t 

then equation 6*2 becomes 

d_P_ _ n P x 
d t x 

Following the procedure in Plow £lj s the above two-* stroke 
standard Diesel cycle processes can be defined in terms of equation 
6*5 as follows. 

1* Compression process* 

He. P e *■ 



fAfe) . _ 

\ d t / rip 



2* Combustion process* 

mi - ° • 



dt / p 

3* Expansion process* 



dtJ 



r\e 



/ 7. ' 

He Pe X 
X 



6.6 



6.7 



6.8 



in the above equation is engine cylinder pressure* 

The electronic analog of the compression or expansion processes 
as represented by equations 6*6 and 6*8 is as follows ■’ 
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1 . 



Integrate (A P e to obtain - P with an appropriate 

initial condition. P • 

P 

2. Multiply -P by n x or n°x to give - n P x or 

7 e J e e a ee 

-nM? x aS a PP ro P r * ate * 

3, Divide by x resulting in the formation of equation 6.6 
or 6.8 and thus returning to the starting point of the 
analog circuit loop. 

The analog of the Diesel cycle applying the above processes and 
procedure is shown schematically in Fig. 6. Also shown is the modifi- 
cation of the circuit of P*low jjLj employing single quadrant electronic 
multiplication for more uniform and stable operation. In this investi- 
gation single quadrant multiplication was employed in each of the cir- 
cuits simulating engine, compressor and bounce cylinder pressures. This 
was accomplished through the use of additional switching relays and 
sign changing amplifiers in order that the electronic multipliers (BEM 
and DEM)* involved in each of these circuits would function only in one 
quadrant . 

A detailed explanation of the circuit of Fig. 6 follows: 

1. The relay switch in the feedback of Amplifier 1 is 

closed whenever x is greater than establishing 

the initial condition, P . 

P 

2. The electronic multiplier and its associated amplifier, 
Amplifier 2, together give a product of one-fiftieth of 
the variable inputs to the multiplier. 

1 

Circuit component symbols defined in Appendix IV, 
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Figure 6* Electronic Analog of the Two-Stroke Standard Diesel Cycle* 






1 






3, Relay 1 is actuated by a change in sign of velocity, 
x* and applies the exponent n^ /2 when the velocity 
is negative and the exponent n V2 when the velocity 
is positive,, 

4* Division of a variable by another variable quantity is 
accomplished by connecting an electronic multiplier in 
the feedback of an amplifier as in the case of Amplifier 4. 
The divisor* -x 9 must always be negative for stable opera- 
tion, while the variable input to Amplifier 4 may have 
either polarity,. The coefficient potentiometer* a^, in 
conjunction with the dividing network* is adjusted to give 
ten times the quotient of the inputs, 

5o Relay 2 is actuated by a change in sign of velocity, x, 
and in conjunction with the relay makes ^ P <g = 0 

d't 

simulating injection of fuel at the commencement of the 
power stroke, 

6, Relay 3, connected in the feedback of Amplifier 6* is also 

actuated by a change in sign of velocity* x* and results 

in Amplifier 6 performing the integration of - 1 n'Px 

10 e 

whenever x is positive. The integration of this quantity 
produces a result which is proportional to the total in- 
stantaneous fuel energy which in turn is compared with the 
desired value of in Amplifier 7„ When the fuel energy 

becomes equal to the desired value as represented by a 
negative voltage corresponding to m^LHV* the relay 
is actuated* reconnecting d P& back to Amplifier 1 

~cHt 
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through Relay 2 0 This action terminatea In* 1 input and 
initiates the polyiio-nic txpan&ion prcc:-'-« 
y o v> 'b^ action of Kul&yo <* mtd 5 is to eftMJfe the use of virile 
qnznrim multiplication in the circuit* IhesMb relays are 
fcuth actuated by a change in sign of velocity, x„ Relay 
4 is so actuated that the input* x, to the electronic 
multiplier is always negative in sign,, Pel&y 5 an 1 Am- 
plifier 5 are therefore required and so connected to ob- 
tain the proper polarity feedback to Amplifier 1 for each 
stroke 0 * 



24 



Elect tonic Analog of the Reciprocating Compressor. 



A typical pressure versus volume Indicator diagram for a recipro- 
cating compressor is represented by r ig c 3(b). The cycle is considered 
to consist of polytropic compression and expansion of air at exponent 
and constant pressure air intake and discharge. 

The computing arrangement for solution of the polytropic equation 
for the compressor, 

d p c = - Pc ■y . 7.1 

dt y 

* 

is similar to that of the engine cycle 0 However, in this case diode 
limiters are employed on the amplifier generating to establish 

limiting voltages corresponding to the intake and discharge pressures. 
This is in conformance with the circuit developed in Plow £l] with the 
exception that single quadrant electronic multiplication was again em- 
ployed involving the necessary two extra switching relays and additional 
sign changing amplifier. 

A schematic diagram of the compressor analog with the single quad- 
rant multiplication feature is shown in Fig. 7. The various functions 
of the circuit components in performing the polytropic compression and 
expansion processes are essentially the same as those described for the 
engine analog In the preceding section. 
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Figure 7. Electronic Analog of the Reciprocating Compressor. 




Electronic Analog of the Bounce Cylinder "Gas Spring". 






Fcr the Bounce cylinder analog, air is alternately expanded and 
compressed by reversible polytropic processes at exponenet result- 

ing in a pressure-volume diagram as shov/n in Fig. 3(d). The differential 
equation involved for the computing loop is 

d Pfr = - H b Pb j , 8.1 

dt 2 

and the schematic diagram of the analog circuit is shown in Fig. 8. 

The employment of single quadrant multiplication in this circuit had 
its most advantageous effect. In the corresponding circuit of Plow 
the absence of this feature resulted in errors which became magnified 
with time by the lack of an initial condition imposed each cycle. This 
modification resulted in stability of the bounce cylinder circuit or the 
ability of the pressure-displacement curve to reproduce or retrace it- 
self over a number of cycles of operation. 
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Figure 8. Electronic Analog of the Bounce Cylinder "Gas Spring’ 



9 e Electronic Analog of the Friction Force. 



In tids investigation a coulomb friction force , "0 n , only was 

considered. This force was assumed to be independent of piston speed 
and gas pressures £7] arid therefore constant in magnitude. An illustra- 
tion of the coulomb friction force is given in Fig. 3(c). 

The friction analog of Plow DJ involving both coulomb and viscous 
friction was modified to produce a voltage corresponding to coulomb 
friction alone and is shown below: 



- X 1 

“VVWW 




Figure 9. Electronic Analog of the Friction Force. 

The limiters in the feedback of the friction analog amplifier 
serve the purpose of establishing the equal magnitude but opposite 
polarity voltages, "C M , corresponding to a constant friction force, 
at the output of the amplifier. With this circuit, employing the 
negative of piston velocity voltage as the input, the appropriate out- 
put polarity voltage for a complete stroke is obtained at the instant 
of commencing the instroke and outstroke of the piston. Those voltages 
are of the proper polarity for the complete analog circuit to be shown 
in the following section. 
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10. Electronic Analog of the Free-Piston Gas Generator. 



The preceding electronic analogs for simulation of the engine, com- 
pressor and bounce cylinder and friction forces were individually assembl- 
ed and tested using the artificial stroke input circuit of Plow flj . 

The artificial stroke consisted of displaced sinusoid voltages simulat- 
ing piston motion at about one cycle in ten seconds based on real time. 
These voltages included the variable M x". M y M , and "z" displacement 
voltages employed in the engine, compressor, and bounce cylinder analogs 
respectively. 

On tying together these component analog circuits with additional 
amplifiers to sum and integrate the various piston forces to obtain its 
motion, the complete free-piston analog was obtained and the artifical 
stroke input no longer required. 

The complete circuit is shown in Fig. 10 and corresponds to the 
basic computer arrangement of Fig. 4. 
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Figure 10. Electronic Analog of the Free-Piston Gas Generator. 



11c Operating Characteristics of the Free-Piston Gas Generator, 

In the free-piston gas generator, the omission of a crank mechan- 

sim permits a free stroke and therefore variable operating conditions 

of piston movement. Since there is no fixed limits for the stroke, the 

motion of pistons should respect certain ranges of inner and outer dead 

point positions, IDP and ODP, as measured from the midpoint of the engine 

cylinder. The limits of ODP should fall within ODP^^ where scavenging 

becomes insufficient and ODP corresponding to mechanical contact at 

the outer limit of piston travel, IDP must lie between the point of 

contact of the two pistons at the midpoint (IDP . = 0) and IDP 

mm max 

where ignition temperature is no longer reached by compression (200 psi). 

IDP . or clearance also should not be so small that excessive peak com- 
mm 

bustion pressures are reached in the power cylinder as compression press- 
ure is dependent on the super-charge pressure and IDP^^. Operating sta- 
bility of a gasifier is thus the ability to keep the piston motion within 
the limits just defined. For the SIGMA GS-34 gas generator the approxi- 
mate operating limits for the ODP are between 16 and 20 inches from the 
midpoint of the engine cylinder. C5] 

Changes in load on a gasifier coupled to a gas turbine are effected 
by adjusting the gas delivery pressure to the power demand on the turbine. 
In general the turbine has no throttle valve and is of a fixed nozzle 
area, full-admission type. The load on a gas generator then corresponds 
to the gas delivery pressure which is related by valve and port pressure 
drops to the scavenge and compressor discharge pressures. 

The frequency of operation is dependent on the gas delivery or load 
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and on the compression pressure in the diesel cylinder. The higher 
these pressures are, the greater is the frequency of oscillation. 

The operating frequency decreases automatically with decreased load, 
or lower supercharge, resulting in lower air flow rate in accordance 
with the decreased turbine requirement. To assist in adapting the gas 
output to changes in turbine load, the position of the ODP is shifted, 
being moved inward for a decrease in power. 

For each pressure level in the gas generator, a certain specific mean 
pressure is required in the bounce cylinders in order to obtain desired 
stroke lengths and compression ratios in the power and compressor cylind- 
ers. For different load conditions, the air pressure in the bounce cylin- 
ders is regulated by a device called the "stabilizer" designed so that 
the average pressure in the bounce cylinders is independent of the stroke 
and varies linearly with scavenge air pressure. Scavenge air pressure 
in turn is dependent on turbine load. The stabilizer therefore provides 
an automatic adjustment of the bounce cylinder pressure, and hence compres- 
sion pressure in the engine cylinder to any imposed operating pressure. 

The compression pressure also may be chosen freely for each delivery 
pressure by proper adjustment of the bounce cylinder pressure. Full 
advantage of free-piston machinery is obtained through use of the stabil- 
izer control which provides for flexible operating conditions. 

In order to vary the gas load of a generator it is therefore necessary 
to change the amount of fuel injected and to change the bounce cylinder 
pressure. Bounce cylinder pressure is dependent on the operating pres- 
sure and is automatically adjusted to new load conditions as previously 
described. 



33 



The power output of a gasifier- turbine combination is controlled 
by changing the delivery pressure of the gasifier and by generating 
gas at the correct rate corresponding to the turbine consumption capa- 
city. However for low power outputs of the turbine, the generator is 
restricted to a minimum gas flow rate which is considerably greater than 
turbine requirements. With reduction of fuel to the gasifier, piston 
frequency and stroke are reduced. This may occur until a minimum stroke 
is reached which must always be sufficient to uncover the ports to pre- 
vent stalling the engine. This operating point defines the "technical 
minimum" of the gasifier and ’according to BARTHALON and HORGEN [5^ 
occurs at 25 per cent of full load for the SIGMA GS-34 gas generator. 
Thus if it is desired to operate a single gasifier-turbine combination 
at a lower output it would be necessary to by-pass the excess gas. This 
of course would result in a direct loss of energy and unfavorable fuel 
economy at low loads. 
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12 0 Simulation of Full and Part Load Performance of the SIGMA GS-34 



Free-Piston Gas Generator on the Analog Computer. 



In this investigation simulation of full and part load performance 
of a gas generator was conducted on the analog computer using the fixed 
geometric and variable operating parameters of the SIGMA GS-34 gasifier. 
The generator was assumed to be directly coupled to a gas turbine with 
an isentropic efficiency of 85 percent. 

Simulation of full load was first performed on the computer with 
results of this run compared with actual performance data [2] and [b] 



followed by part load runs of 75, 50, and 25 percent of full load (25 
percent load being so called M technical minimum") 0 Gas delivery pres- 
sures, P^ , for these loads were obtained from London £,2J and Barthalon 
00 and are as follows (refer to section 5 for system parameters): 

Load Percent Gas Delivery Pressure, P^, (psia) 



100 


64.5 


75 


54.0 


50 


45.0 


25 


37.0 



As related in the previous section, mean effective bounce cylinder 
pressure was considered to vary linearly with scavenge air or gas deliv- 
ery pressure. Eichelberg [3j also shows essentially a constant pressure 
differential between maximum bounce cylinder pressure and scavenge air 
pressure over the entire operating load range of the SIGMA GS-34 gas 
generator. Thus with initial or maximum bounce cylinder pressure avail- 
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able for full load, [4] values of this pressure may be determined for 
any operating condition* 

The analog computer circuit established in this investigation simu- 
lates a steady state operating condition of a gas generator and is in- 
capable of solving a transient problem such as start up of the gasifier 
or a change in load condition. The initial condition for the computer 
problem is with the piston at the outer limit of its stroke (ODP) or 
initial piston displacement (x q ) and the bounce cylinder pressure at its 
maximum or initial value (p£). A cycle of operation would involve the 
travel of the piston inward to’ the inner dead point position (IDP) and 
return to the initial starting point. Section 3 describes in detail the 
dynamics and forces involved in the piston motion. 

With the computer properly programmed for simulation of a particular 
load as set down in detail in Appendix II the criteria for performing 
each run was to vary initial piston displacement position, x^, so that a 
maximum engine compression pressure of close to 1500 psia was attained at 
the end of the instroke and to adjust the amount of fuel input to return 
the stroke to the initial starting point. The figure 1500 psia was select- 
ed as a nominal value for the standard Diesel cycle employed in this in- 
vestigation considering normal operating values of compression pressure 
for the SIGMA GS-34 gas generator. [3] The above adjustments of initial 
piston displacement and fuel setting had to be made while observing 
engine pressure versus piston displacement diagram on the recording equip- 
ment. With these adjustments properly made, two to three fairly reproduc- 
ible cycles could be attained representing the steady state operating con- 
dition for the particular load run. The adjustments at the same time were 
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rather critical due to inherent drifts in the computer circuit compon- 
ents. 

Comparison of analog computer solution results with actual per- 
formance data for the SIGMA GS-34 gas generator at full load condition 



is as follows 0 . 




GS-34 


Analog 


1. Gas Turbine 


Power output, shp 


1138 


1290 


Gas flow rate, lbm/hr 


29000 


39400 


Gas pressure at inlet, psia 


64.5 


64. 5 1 


Gas Temp, at inlet, 


deg F abs 


1405 


1174 


2. Gas Generator 


Frequency, cpm 


613 


649 


Stroke, inches 


17.5 


18.25 


Engine effective stroke, inches 


9.7 


10.0 


Engine clearance, inches 


1.29 


1.0 


Compressor clearance, inches 


2.38 


2.0 


Bounce clearance, inches 


8.0 


7.7 


Compressor air delivered, 


Ibm/ cylinder- cycle 


0.388 


0.507 


Ratio engine/compressor air 


0.47 


0.404 


Air- fuel ratio 


32 


35 



^Inputs to the computing problem included actual performance data 
pressures as here shown in addition to fixed geometric parameters. 
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GS- 34 



Analog 



Engine Indicated thermal 

efficiency, LHV basis 29.8 51.3 

3. Gas Generator-Turbine System 

Fuel flow rate, lbm/hr 460 452 

Brake specific fuel consumption, 
lbm/shphr 0.404 0.351 

Thermal efficiency, LHV basis 34.6 39.8 

On a net work per cycle basis, engine work is equal to the sum 
of compressor and friction work from equation 3.3 and for the computer 
results at full load a discrepancy of 17.0 percent on basis of net 
engine work was obtained. 

Figs. 11, 12, and 13 show the full load pressure versus displace- 
ment diagrams for the engine, compressor, and bounce cylinders re- 
spectively. 
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Pressure* psia 




Figure 1 1 « Pressure versus Displacement Diagram for 
Engine Cylinder at Full Load. 



39 



Pressure f psia 






100 
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Figure 12. Pressure versus f.isp lacement Diagram foi Compressor 
Cylinder nt Full toad. 




JC 25 20 15 10 5 

Displacement , inches 

Figure 1 . 1 re ; ->ure versus Displacement Diagram for Bounce 

Cylinder at Full Load. 
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Pressure, psia 



Computer solution results for part- load performance with operating 



pressure varied as previously described are shown in the table below 
for load percentages as indicated: 



lo Gas Turbine 

Power output, shp 
Load percent on basis of 
computer full load results 
Gas flow rate, lbm/hr 
Gas pressure at inlet, psia 
Gas temp* at inlet, deg F abs 
2, Gas Generator 

Frequency, cpm 
Stroke, inches 

Engine effective stroke, in. 
Engine clearance, in. 
Compressor clearance, in. 
Compressor air delivered, 
lbm/cyl-cycle 
Ratio engine/comp, air 
Air- fuel ratio 
Engine ind, thermal eff., 

LHV basis 



Percent 


of Full 


Load 


75 


50 


25 


963 


687 


442 


74.7 


53.2 


34.2 


36100 


33800 


31000 


54.0 


45.0 


37.0 


1066 


925 


767 


605 


561 


517 


18.05 


17.85 


17.65 


10.1 


10.15 


10.15 


0.9 


0.85 


0.85 


2.0 


1.95 


1.95 


0.497 


0.503 


0.500 


0.360 


0.309 


0.268 


36.7 


40.5 


51.0 


0.516 


0.572 


0.656 
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Percent of Full Load 



3. Gas Gen.-Turb. System 


75 50 25 


Fuel flow rate, Ibm/hr 


354 258 163 


BSFC, lbm/shphr 


0.368 0.375 0.369 


Thermal eff., LHV basis 


38.0 37.2 34.3 



Net work discrepancies between engine and sum of compressor and 
friction work for the part-load conditions are as follows: 



Percent of Full Lbad 


Discrepancy, percent 


75 


9.0 


50 


4.6 


25 


7.2 



Figures 14 throxigh 22 show pressure versus displacement diagrams 
of the engine, compressor, and bounce cylinders for the three part- 
load conditions* Figures 23 through 29 show performance characteris- 
tics for full and part-lead operating conditions as obtained from the 
analog computer solution results* Also 9 shown on these graphs are the 
actual performance operating points for the GS-34 gasifier for full 
load condition* 
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Figure 14. Pressure versus Displacement Diagram for 
Engine Cylinder at Three-Quarters I,oad„ 
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Pressure, psia 




Figure 15. Pressure versus Displacement Diagram tor Compressor 
Cylinder at Three-Quarters Load. 
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Figure 16. Pressure versus Displacement Diagram for Bounce 
Cylinder at Three-Quarters Load. 
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Pressure, psia 



Pressure* psia 




Figure 17c Pressure versus Displacement Diagram for 
Engine Cylinder at One-Half Load. 
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Figure 18. Pressure versus Displacement Diagram for Compressor 
Cylinder at One-Half Load. 
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Figure 19. Pressure versus Displacement Diagram for Bounce 
Cylinder at One-Half Load. 
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Pressure* psia 



Pressure* psls 







Figure 20. Pressure versus Displacement Diagram for 
Engine Cylinder at One-Quarter Load. 
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Pressure, psia 




Figure 21. Pressure versus Displacement Diagram for Compressor 
Cylinder at One-Quarter Load. 
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Figure 22. Pressure versus Displacement Diagram for Bounce 
Cylinder at One-Quarter Load. 
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Pressure, psia 



Thermal Efficiency, percent 




Figure 23. Thermal Efficiency on LHV Basis of Combined Gasifier- 
Turbine System versus Load from Computer Solution 
Results. 
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BSFC, lbm/shphr 
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Figure 24. Fuel Consumption versus Load from Computer Solution 
Results. 



50 



Fuel Rate, lbm/hr 



Turbine Inlet Pressure* psia 
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Figure 25. Turbine Inlet Pressure and Temperature versus Load 
from Computer Solution Results. 
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Turbine Inlet Temperature, °F abs 



Frequency, cpm 
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Figure 26. Frequency and Gas Flow Rate versus Load from Computer 
Solution Results. 



52 



Gas Flow Rate, lbm/hr 
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27. Operating Pressures versus Load Employed in Computer 
Simulation of SIGMA GS-34 Gas Generator. 



IDP, inches 




Figure 20 o Inner and Outer Dead Points as Measured from Midpoint 
of Engine Cylinder versus Load from Computer Solution 
Result So 
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Figure 29. Ratio of Air/Fuel and Engine/Compressor Air versus 
Load from Computer Solution Results. 
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Air/Fuel Ratio 
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13* Conclusions. 



The employment of single quadrant multiplication in the various 
computing loops resulted in uniform and stable operation of the analog 
circuit as a whole* However* the lack of automatic balancing permitted 
a certain amount of amplifier drift which could account for appreciable 
errors during the time required to complete a particular load run. 

Engine, compressor and bounce cylinder indicator diagrams had to be 
taken individually on the single available X - Y variable plotter. 

In comparing computer solution results with actual performance 

* 

data for the SIGMA GS-34 gasifier at full load condition (section 12), 
consideration should be made of the various assumptions and idealizations 
in constructing the analog circuit, particularly in the type of engine 
cycle employed. 

Quite possibly a further refinement in this type of an investiga- 
tion would be to adjust appropriate circuit parameters such as poly- 
tropic exponents to make the full load analog results come in closer 
agreement with actual performance data. With the resulting analog 
more closely simulating a given gas generators operating characteristics, 
a more accurate analysis or prediction of part load behavior could be 
expected. 

The performance data obtained from computer solution results for 
the various load conditions as tabulated in section 12 and graphically 
represented in Figs. 23 through 29 are considered reasonable for the 
operating conditions or criteria established for making the runs. The 
flexibility of the free piston system would permit a large variation of 
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running conditions for a given output . 

In the full and part load performance studies of the gas genera- 
tor-turbine system of the previous section* no consideration was made 
of the flow characteristics of the turbine coupled to the gasifier. In 
actual practice to avoid losses by throttling between the gasifier and 
turbine it would be necessary to match the delivery of the gasifier 
with the consumption of the turbine. This is required since, with a 
given turbine, for each gas pressure there corresponds a particular mass 
flow. A particular mass flow requirement at a given delivery pressure 
from a gas generator would fix its stroke length. 

This investigation has demonstrated that the operation of a free 
piston gas generator may be simulated to a reasonable degree of accuracy 
and stability by electronic analog means t By more closely approximating 
the various cycles of the actual machine, it is considered that the 
analog method could be an important means of predicting the performance 
of free-piston engine systems. 
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APPENDIX I 



TIME AND MAGNITUDE SCALING 



With the exception of the scaling for the friction analog, the 
same time and magnitude scaling relations as developed in Plow £1J apply 
to the computer circuitry of this investigation. These scaling relations 
with respect to the circuit components of Fig. 10 are given below: 

1. Engine pressure, compressor pressure, bounce pressure and 
displacement integrating amplifiers (Nos. 1, 8, 12 and 19 respectively). 






Asl 



= <3l<? = ^ X 



Within the computing loops simulating engine* compressor, and bounce 
pressures, magnitude scales involved in the multiplications and divisions 
are fixed as described in sections 6, 7 and 8. Scaling is required only 
for the above pressure integrating amplifiers. 

2, Summation of forces amplifier (No. 17) 



and 





Ae °<Pe . 


®im 


M * x 


&11& %-fn - 


Ac Pc . 


&11Z 


M °<-y 


3 R-Piq = 




^193 


M «x 




•^F + r . 
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3. Velocity integrating amplifier (No. 18), 

&18 = X . 

t~fi8 ^x tX'-b 
60 



The above scaling relations fog&th^r with tlxvd circuit qurintl • 
ties as shown in Pigs. 6 9 7* 6 and 9 aie all that ar* resulted ior 
thtt complete analog network of Fir;. 10 o 



ATl-siri'ir:- T. 

I' Kl\, RAMMING AGALvo r " HT ..il! G £ * ' '. + 

FREE-PISTON GAS GENERATOR - FULL AtJi P.v a 

Programming of the analog computer for the STGKA GG-iA free- 
piston gas generator for full and part load conditions requires the 
following input parameters: 

1. Fixed geometric parameters (references 2 and 4): 

0.979 ft 2 
5.73 ft 2 
6.77 ft 2 
1110 lbm 
11 .0 in 
1.09 in 
26.92 in 

2. Fixed operating parameter (section 3 and reference 2): 

C 960 Ibf 

3. Variable operating parameters (references 2 S 4 and 6 and 
section 13): 

Percent of Full Load 

100 75 50 25 



p i’ 


psia 


14.0 


14.0 


14.0 


14.0 


V 


psia 


67.9 


56.7 


47.2 


38.9 


V 


psia 


64.5 


54.0 


45.0 


37.0 


p b> 


psia 


132.3 


121.8 


112.8 


104.8 



A 

e 



c 

*b 

M 



x 

P 

c 



b 
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4« Assumed polytropic exponents 0 



n 


1.30 


c 


n' 


1.40 


e 


n b 


1.40 


ri 


1.30 



c 



The various scale factors required for the problem are assigned 
values as follows^ 



p £ ss 50 psia/volt * 


7200 psia/volt 


oc P( . a 5 psia/volt = 


720 ps fa /volt 


P =£ 8 psia/volt s 


1152 ps fa/volt 


^ Ffr~ lb£/volt 




©4.^ ~ 0.5 in/volt - 


1/24 ft /volt 


o< ^ - 2 ft/sec/volt 




c< * s 100 ft/secVvolfc 




^ = 100 sec computer 


time/i sec tea 



Based on the above scale factors^ computer voltages corresponding 



to the geometric and operating input parameters and ss has: red 



quantities are listed belows 

x 22 o 0 

P 

c 2« 18 



b 

C 



P 

P 

T 

i 



1 



d 



P 



o 

b 



53o 84 
12 0 0 



Percent 


of Full 


Load 




100 


75 


50 


25 


1.29 


1.08 


0 . 90 


o.74 


2.8 


2.8 


2.8 


& 


13.58 


11.33 


9,45 


7.77 


16.55 


15.25 


U.l 


l ".1 


63 









Circuit components .".re determined from the scaling relationships 



of Appendix II as follows: 

1. For pressure and displacement integrating amplifiers 

(Nos. 1, 8, 12 and 19), 

3l * S_8 - = 3 iq » = 0.480 

R l C fl ^8 Cffi ^12. Cfj.2. Riq C-fiq °<X 

R's = 1 meg 

C's - l/<fd 

a’s 0.480 

2. For summation fof forces amplifier (No. 17), 

3m n s 4g ^Pe. a 2.044 

fcini M •<■ * 

<^172, = Ac ** • Pc = 1.196 

R 1 7 2. M °c x 

3 173 R-fl7 « = 2.262 

Rl73 M°<X 

3p4 + r = 0.0232 

R 1 7 4 M 04 x 



f 17 


= 


1 meg 






171 


S 


0o25 meg; 


*171 


« 0.511 


R 172 


S 


0o 5 meg; 


a 172 


* 0.598 


173 


SS 


0*25 meg; 


a l73 


- 0.565 


R 174 


= 


10 meg; 


d 174 


- 0.232 
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3 0 For velocity integrating amplifier ** 18) 

m X C 0, 560 

&1& C-f 18 c *’b 

R lg ~ 1 rr<e 8 
C f 18 “ 

a 15 = 0.500 
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APPENDIX III 



FULL AND PART LOAD PERFORMANCE CALCULATIONS 
FROM COMPUTER SOLUTION RESULTS 

Data from the computer solution results as obtained from the 
recording equipment for the various load runs is as follows: 

1. For full load 

From time recorder trace 

Cycle period, computer time * 9.25 sec/cycle 

From engine cylinder diagram. Fig. 11 
« 38.5 volts 

X* 3 a 2.0 volts, P 2 = 32.0 volts 

x* 4 ■ 4.5 volts 

P,j ■ 3.7 volts 

2 

Area ■> 144.0 volts 

Prom compressor cylinder diagram, Fig. 12 
y*j « 41,0 volts 

y"^ m 4.0 volts 
■y^ ■ 11.2 volte 

2 

Area «■ 198,4 volts 

From bounce cylinder diagram, Fig. 13 

* 15.4 volts, Pj a 132,3 volts 

i *2 * 52,0 volte, P^ « 3.5 volts 

P* * 7.13 volts 

ave 
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2. For three-quarters load 



From time recorder trace 

Cycle period, computer time » 9.9 sec/cycle 

From engine cylinder diagram. Fig. 14 
= 38.0 volts 

x^ - 1.8 volts, P 2 = 32.0 volts 

x 4 a 3.9 volts 
P^ - 3.0 volts 

Area s 121.6 volts 2 

■a 

From compressor cylinder diagram, Fig. 15 

y*. » 40.2 volts 

5T - 4.0 volts 

F 4 - 11.0 volts 

Area - 171.2 volts 2 

From bounce cylinder diagram. Fig. 16 

T x - 15.8 volts, P^ s 15.25 volts 

z~2 ~ 52.0 volts, P 2 = 3.2 volts 

P *s 6.64 volts 
ave 

3. For one-half load 

From time recorder trace 

Cycle period, computer time =■ 10.68 sec/cycle 

From engine cylinder diagram. Fig. 17 
x = 37.4 volts 

5T = 1.7 volts, P" 2 - 33.0 volts 

- 3.3 volts 

4 

P - 2.2 volts 

5 

2 

Area * 105.6 volts 
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% 








From compressor cylinder diagram, Fig. 18 

y"j = 39.6 voles 

y" 3 - 3.9 volts 

= 10.0 volts 

Area = 155.2 volts 2 

From bounce cylinder diagram. Fig. 19 

T = 16.0 volts, as 14.1 volts 

Z 2 - 52,0 volts, - 3.0 volts 

F =6.13 volts 
ave 

4. For one-quarter load 

From time recorder trace 

Cycle period, computer time ~ 11.6 sec/cycle 

From engine cylinder diagram. Fig. 20 
= 37.0 volts 

= 1.7 volts, P 2 = 30.2 volts 

- 2.9 volts 
4 

P^ = 1.5 volts 

2 

Area = 83.2 volts 

From compressor cylinder diagram. Fig. 21 
Fj = 39.2 volts 

jT-j = 3.9 volts 

y^ = 9.8 volts 

2 

Area = 115.2 volts 

From bounce cylinder diagram. Fig. 22 

- 16.2 volts, Fj = 13.1 volts 

F - 52.0 volts, - 2.8 volts 

P ^ 5.64 volts 

ava 
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Calculation of performance teaultt. f.to<« tail i > »d ate 



foil*, vr 



(Xj - x 3 )c< x * (36.5 



2 . 0 ) 0. 5 



TJo< 

.3 x 



1 . Frequency 

Cycle period based on computet time 

- 9.25 sec/cycle 
Frequency based on computer time 

» 6.49 cycles/mio 

Frequency based on real t line 

- 6.49 X 100 - 649 cycles/min 

2. Stroke and clearances 
Stroke 

- 18.25 inches 

Engine effective stroke - (x^ 

= ( 22.0 

= 10.0 inches 

Engine clearance = x. ' ’oc ) - 2.0 X 0.5 

X 

» !o0 inches 

Compressor clearance 

= 2 o 0 inches 

Bounce clearance « z, > ^( c K f ) 15*4 X 0 o 5 

1 z 

i 

^ ? 0 7 inches 

3o Compressor discharge air 
From equation 5«4 0 4s 

Pi A c Yi 



2 . 0 ) 0. 5 



y 3 <°q> 



5 X C.5 



mJL = 



R. To 
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pi = (4,0 pSi(X 

A c = 5"* 7 3 -ft 2 

Y i - ( y i ~ O4) 

- (4/.0 - 11 . Z) o.S~ s 14 . q in. 

T 0 - 5T5 0°R 

0 ± ■QY/44X5',73y/4. < ?) 

(!T3.3XS‘30X/Z) 

w 0,5"07 /6m 



4. Engine intake air 



From equation 5.4,2: 




Pp - 6? 4 . 5" p S l O, 

w _ (11* OX 0411) 
Vp - 

^7.9 pstfl. 

Pi. * 14*0 psio. 



0,897 4^ 



n c = /. 3o 

m = U+>sXi44Xo-S 12 I 

e ( 6“3 « 3 X 5*3 o)(4 , 8 5"J °* 3< M.3 o 

= 0, 2.0 5T /6m 



5. Ratio of engine to compressor air 

m /m, a 0.205/0.507 
e d 

= 0.404 



6. Fuel input 

Q-f = ^+ LHN/ = St P 3 (V 4 -V 3 ) 

P 5 * P' B (°<Pe) = 3 Z .0 k S’O 

~ \<oOO pSL (k. 



•'O 



A V = 



fo.- X^°<x)Me) - (A.s-z.oYo.rXo.li'i) 

12. ~ 1Z. 

= 0,102. -ft 3 

<?* = m f L HV « (lb0oXH4Xo.loZ) 

= 82., ZOO -ft - Ib’f 

Assuming LHV = 18,200 ETU/lbm 

mx - 8^/00 

+ (i8 / z.ooX r ne) 

*■ 0 , 005‘8 ibny'cycle. - s i de. 

7o Air - fuel ratio 
•% 

A/F = m /m . = 0. 205/0,0058 
e f 

= 35.3 

8. Fuel flow rate 

Fuel rate = 2(m^.) (fr) (60) 

= 2(0.0058) (649) (60) 

= 452 lbm/hr 



9. Gas flow rate 

Gas rate = 2(m ) (fr) (60) 

d 

= 2 (0.507) (649) (60) 

= 39,400 lbm/hr 



10. Engine exhaust gas temperature 



From equation 5.5.1; 

Te .i[l + (ft-l)£] 

T _ _ p s~ v*- _ (i8r)(J4tx o.eqn) 

& ^cR ’ Lo,Z.05)(S$,$) 

£, 1 8 5~ 

1.40 



“TV - 



£i+ (1.40 



= 1 , 1 1 8 °R 



Z t J8S°R. 



n 



lie Generator gas delivery temperature 



T-t - 



From equation 5.5.3: 

T e 4- (l - ' Vle - 



)? 



PA \ ft r — 1 ^ , 0,3 0 / 

' * 5~3o(4.8tr) ^ l0 



T-t - 



Tp = T 0 (|f) 

= U3 0 £. 

0,404 (l, 118)4 (1-0, 404) 103 

1)114 °K 



12. Power output of the system 
•% 

From the gas tables - air at low pressure for one 

pound - isentropic available energy, h , on expansion from state in 

s 

item 11 to atmospheric pressure is equal to 97.8 BTU/lbm. 

From equation 5.5.5: 



. (^^iXT r >8)(Z)(mdX-frXKt) 

(Shp 'turbine “ 33,000 

( c il'6)(in8)(Z)(0.5-07)(6>4q)(0,S5-) 

33,000 

= 1290 hp 

13. Net work balance per engine cycle 
From equation 3.3: 

Wfe ne t erl g. = Wfe „e-£ com p +■ ^^-fr 

w *-neten$. ' 044 -oy^zooXo, 914X^4) 

- 41, zoo -f'b-lb'f 

w kne^comp = (M 8 ,4Yl ZoXS.ISX-Jt) 

= 34 t ioo -ft- Ib-f 

Wkir * ( Z-Xl <= OXS8 , Z$rfz) 

= Z^ZO -fi- lb-f 
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Discrepancy: 7180 ft-lbf or 17.0 percent on 



basis of net engine work. 

14. Brake specific fuel consumption 

BSFC = m / (shp) 

f v ^'turbine 

= 452/1,290 

= 0.351 lbm/shphr 

15. Engine indicated thermal efficiency, LHV basis 

/ w \ t ^ ^ r\€'i) en q i n e .. ^ 

(^engine" 



4£ f ZOO 
82, ZOO 



x loo 



- SI . 3 percent 

16. Overall thermal efficiency, LHV basis 

/ S (2hpX3±£°oXloo) 

V V-t^/ove rail “ 



("VtoX'-^Xvis) 

_ (1,2^0)03,000X^00) 
(tsz/ 60 )(lS ) 2oo)(nn8) 

— 3*3 .8 percent 



Above sample calculations represent the full load condition of 
the gas generator- turbine system. Performance results for part- load 
runs are tabulated in section 12 of the investigation. 
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APPENDIX IV 



ANALOG COMPUTER SYMBOLS 



Circuit 

1. Potentiometer: 



ez 



Symbol 



-o- 

a 



Function 



So = 36 



2. Summer; 



_ ft. 

I-WW-, 



o. Rz 

C ^<>A/WV- 



a 3 R = 




6l Q^l 

eiod2j R i 

^3 q 33 




_^_4 



!4=-(^e 1 + ^e ; 

.) 



ei + 

*3 



Summing Integrator: 




e 4 =- 









e, + 



+ ^-e 3 



r 3 c^> 



c-f^ ' RzC + ei 

,j JL-b + E c 



4. Relay Driver: 




0 w Ae n 4. 3 a. 6 7. n, 

Rx 

(+6 uohevhJLi + atet < 0 

R i R-i 
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Circuit 



Symbol 



Funct ion 



5* Boeing Electronic Multiplier: 



10 







-L- p p 

SO "1 C 2. 



6 0 Donner Electronic Multiplier: 








7. Division using Boeing Multiplier: 







— 


e«. 








BEM 








BEM 















^ — O'W/v 




-e*. 







e 3 






- so 



di£i 



8. Division using Donner Multiplier: 




In division circuits using Boeing and Donner mult if-1 ters s variable 
divisor voltage, e^, must always be negative for stable operation. 



75 




APPENDIX V 



DESCRIPTION OF EQUIPMENT 

For this investigation the following equipment components were 
employed: 

1. Two Boeing electronic computer (BEAC) racks D-\) contain- 

ings 

a. Operational amplifiers . 

b * Diode limiters* 

c* Electronic function multipliers* 

d, Coefficient potentiometers* 

e, I* C* voltage source* 

2* Donner electronic function multipliers* 

3* Sanborn time recorder* 

4* Electronics Associates X-Y plotter* 

5* Electro Instruments digital voltmeter* 

6. Single pole, double throw, non-polarized sensitive relays* 

The above equipment is shown in Fig* 31 . 

The computing network for simulation of the free-piston gas genera- 
tor problem required a total of 25 operational amplifiers, 6 electronic 
function multipliers, 11 diode tube limiters and 10 sensitive switching 
relays* 

The switching relays were actuated by the output of an operational 
amplifier* Amplifier feed-back arrangement and connection of the relay 
energizing coil are shown in Fig* 20 next page* 



76 




Figure 30. Actuation of Relays. 

Sensitive relays available for this investigation were of the 
single pole, double throw, non-polarized type. With no energizing 
voltage, one of the two sets of contact points would be closed by 
spring action, i/hen actuated, contact between the other set of points 
would be made against the spring force. 

As indicated in Fig. 30 the feedback of the relay driver amplifier 
consisted of diode limiters and a bias voltage, £ . A positive voltage 
signal at the grid of the amplifier would effectively clamp tine output 
of the amplifier to ground potential resulting in an untctuated condition 
of the relay. A negative input of any magnitude on the other hand would 
result in an amplifier output determined by the bias voltage thereby 
actuating the relay. Bias voltage was adjusted to a value sufficient 
to operate the relays but not to exceed the maximum 10 milliampete rat- 
ings of the amplifiers. 
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% 














The above principle was employed tor the actuation of all relays 
identified by the subscript "R" in the computing circuit of Fig. 10. 
In this figure the eight x relays (energizing coils all in series) 
were actuated by Amplifier 21, x p relay by Amplifier 23 and Q f relay 
by Amplifier 7„ 



78 



Electronic Multipli 
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